Ratchet effect in dc SQUIDs 
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We analyzed voltage rectification for dc SQUIDs biased with ac current with zero mean value. 
We demonstrate that the reflection symmetry in the 2-dimensional SQUID potential is broken 
by an applied flux and with appropriate asymmetries in the dc SQUID. Depending on the type 
of asymmetry, we obtain a rocking or a simultaneously rocking and flashing ratchet, the latter 
showing multiple sign reversals in the mean voltage with increasing amplitude of the ac current. 

' Our experimental results are in agreement with numerical solutions of the Langevin equations for 

, the asymmetric dc SQUID. 

O 
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Directed molecular motion in the absence of a directed net driving force or a temperature gradient in biological 
, systems has drawn much attention to the so-called Brownian motors |^J^ where nonequilibrium fluctuations can 
• induce a net drift of particles along periodic structures which lack reflection symmetry. Such structures are denoted 
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as ratchets. 

In general, a ratchet system is characterized by the equation of motion of an overdamped particle in a periodic 
\0 ' potential 

o ■ 

O' ^^ = -^[w{x,t)-Fd{t)x]+F^{t)^-^u{x,t)+F,,{t) (1) 

o : 

■ with w{x,t) ^ w{—x,t) for all choices of origin. Here, ^ is a friction coefficient, x a cyclic coordinate, Fd{t) a driving 
force and Fn(t) represents a force with zero mean value due to thermal noise with Gaussian distribution. One can 
classify different types of ratchet systems depending on the actual form of w{x, t) and Fd{t) (see e.g. 1^). The rocked 
thermal ratchet |^|^ is obtained if one chooses w{x,t) to be time independent and Fd{t) to be either stochastic or 
deterministic. For the flashing ratchet a time dependent potential w{x, t) is chosen which, in the simplest case, is 
^ ■ of the form w{x,t) = wo{x)wd(t), where Wd{t) again can be stochastic or deterministic. The common feature of all 
Q ' ratchet systems is their rectifying property: it is possible to extract directed motion from nonequilibrium fluctuations 
O ■ or periodic excitations with vanishing time average. Although considerable theoretical work exists in particular for 
^ ^ stochastically driven ratchets, only few experiments have been realized e.g. by means of dielectric and optical 
potentials 0. 

The periodicity of the Josephson coupling energy Ej with respect to the phase difference S of the macroscopic 
wave function across a Josephson junction is an ideal prerequisite to build a ratchet system. Within the resistively 
and capacitively shunted junction model the equation of motion of 6 is equivalent to the motion of a particle in 
the so-called tilted washboard potential U{6,t) = —^{IqCOsS + I{t)S} (see e.g. |^). Here, $o is the flux quantum 
and Iq is the maximum Josephson current across the junction. An external force causing the tilt is obtained by a 
bias current /, and the mass and friction coefficient are determined by the junction capacitance C and the junction 
resistance R, respectively. In the strongly overdamped limit {Pc = ■^IqR^C ^ 1) C can be neglected. The voltage 

V — across the junction represents the velocity of the particle via the second Josephson equation. For a short 

Josephson junction with a current phase relation which is an odd function, U{S) is always an even function and no 
ratchet potential is obtained. However, voltage rectification may be induced by asymmetric fluctuations as proposed 
in 1^]. To obtain a ratchet potential for Josephson junction systems, two ratchet configurations based on coupled 
Josephson junctions have been proposed. Zapata et al. | |lO| considered an asymmetric three junction dc SQUID, with 
vanishing loop inductance L thereby coupling the phase differences across the junctions rigidly via the applied flux 
<I>Q. Due to the doubled phase shift in the SQUID arm with the two junctions in series the reflection symmetry in 
the effective 1-dimensional potential is broken. Falo et al. proposed a 1-dimensional Josephson junction array 
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with spatially alternating critical junction currents and alternating loop inductances. They showed that a Josephson 
fluxon experiences a ratchet potential as it moves along the array. This effect has been confirmed experimentally very 
recently |l2). 

In this letter we show that the well known voltage rectification in asymmetric dc SQUIDs jl^ is a consequence of 
broken symmetry in the two dimensional SQUID potential due to this asymmetries and the application of magnetic 
flux. In contrast to the model proposed by Zapata et al. this results in a quasi two dimensional ratchet potential. 
We demonstrate its experimental realization by using a two junction high transition temperature (Tc) dc SQUID 
with small but finite L. In our experiments we use an asymmetry aj in the critical currents /q — (1 — a/)/o and 
/q = (l + a/)/o of the two junctions and/or an asymmetry ul which introduces a self-field effect, i.e. afiux component 
'^seif = ocl'^qI / Ic with /c = /g + /q. a finite value of aL can be realized for example with different inductances in 
the left and the right arm of the SQUID loop; our definition however is more general and accounts for all kinds of self 
field generation. 

The extension of the tilted washboard potential model in two dimensions has been used to provide insight into the 
dynamics of both, symmetric dc SQUIDs and dc SQUIDs with an /q asymmetry ||l5). The SQUID potential Us 

Us{5u52,t) = -^{ll cos + /2 COS 62 I . 

-l^iS^ - ^1 - + + + 62)} 

contains three basic terms, (i) the two cosine terms representing Ej, (ii) the quadratic magnetic energy term and (iii) 
the external force term proportional to /. Here, 61,62 are the phase differences across the two Josephson junctions 
and /3l = IcL/^o is the normalized inductance. An asymmetry aj ^ introduces a difference in Ej for the two 
junctions thus deforming the double cosine part of the potential landscape, while aL 7^ shifts the minimum of the 
parabolic magnetic energy term proportional to /. Transforming into more appropriate phase coordinates along the 
direction of the bias current 6 = {61 + 62)/'^ and the applied flux ip = {62 — 61) /2 and using i = I/Ic and fa = ^ 
yields in analogy to eq. (|^) Us{6, ip, t) — 'Ws{6, ^p, t) — i{t)6 with 



'Ws{6,ip,t)= — cos (5 cos (yj + a/ sin (5 sin 
1 



+^{^-■^/a-7^ai^(^)}2. (3) 



In this coordinates, the voltage V across the SQUID is determined by A ratchet effect, i.e. voltage rectification 
for an external excitation with zero mean value, will occur if the system is not invariant under parity transformations 
with respect to 6, that is for 

Ws{6o + 6,ipo + ip) ^ Ws{6o - 6,(pQ + ip) and 
Ws{6o + 6,'po + p) ^ Ws{6o ~ 6,'po - (p), 

for arbitrary origin {60, (po)- Eq.(§) satisfies these conditions, i.e. Ws is a ratchet potential in (5, if /a ^ f A {aj ^ 
V aL 7^ 0); (n = 0, 1, 2...). For given a/, a^, the applied fiux fa is an experimentally controllable parameter, which 
allows to change the parity of Ws. The type of asymmetry of the dc SQUID driven by a fiuctuating or periodic bias 
current i{t) determines the nature of our ratchet system. For aj 7^ 0, = 0, Ws is time independent, thus we have 
a rocking ratchet system. For aL ^ 0, Ws becomes time dependent since the minimum of the magnetic energy term 
shifts proportional to i{t). In that case we obtain a simultaneously flashing and rocking ratchet. 

Figure l] shows contour plots of the 2-dimensional SQUID potential Us for = and different values of ai and 
bias current i. For ai — (upper row), the potential has reflection symmetry with respect to 6, and with finite i we 
have Us{6, ip; i) = Us{—6, ip; —i) and thus V{i) = —V{—i). For aj — 0.3 (lower row) no reflection symmetry is present 
and in fact, while for i = +0.9 no local minima are found i.e. the solution for the phase is unbounded, for i — —0.9 
local minima are present and no voltage appears. 

The signature of a ratchet effect can be observed in the V{^a) characteristics of an asymmetric SQUID (a/ ^ 0) 
in form of a fiux shift A<I> (see Fig. ^) between the two curves obtained by biasing the SQUID with the currents ±/ 
with an appropriate choice of |/| « Ic (experimental details are given below). A flux shift ^ n$o corresponds to 
V(l, '^a) 7^ — ^(— -f, '^a) for all 7^ ^$0 ("■ =0,1,2..). Thus the bias current (or time-) reversal symmetry is broken, 
which is true for every tilted ratchet potential: for certain values of the external force the magnitude of the velocity 
of the system depends on the sign of the force. Note however, that y(J, <I>q) = —V{—I, — $a) always holds. Applying 
an ac excitation at angular frequency uj with zero mean to a SQUID with A<I> 7^ n<I>o results in a finite mean voltage, 
which is maximum for <I>a = <^i^ and zero for <I>a = §$o- Obviously the maximum mean voltage with ac bias is 
largest if A$ = 0.5$o- 

Next we analyze the dependence of A$ on a^, a/, /3l and /. With self field effect and symmetric junctions (cth 7^ 0, 
cti ~ 0), A<i> between a pair of V^($a; ±^) curves is simply proportional to the bias current and a^'- A$ — ^^^a^l j Ic- 
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For aL = 0, a/ 7^ 0, A$ depends on a/, Pl and can be obtained as follows. The flux value <I>i (cf. Fig. |[) for / Ic 
is equal to the flux value of maximum critical current Ic = /q + /q . Ic is reached if the circulating current in the loop 
equals a/Jo, corresponding to an applied flux ^^oc>:iPl- From our definition it follows that A$ — $oQ!//3l- For a 
SQUID with given inductance this relation can be used to determine the junction asymmetry ai by measuring A<i> 
and Ic- Wc note that A<i> 7^ independent on / is frequently observed for our high-Tc SQUIDs which is not surprising, 
since the spread of the critical currents of the two junctions can be substantial ||l6|] . 

To obtain more insights into the dynamics of the system we performed numerical simulations to solve the coupled 
Langevin equations of the asymmetric dc SQUID [l^ with Pc ^ 1 and with thermal fluctuations due to Nyquist noise 
in the junction resistors which adds white noise with zero mean and Gaussian distribution. We use the normalized 
noise parameter F = '^J}^^ with thermal energy fc^T. In our simulations we chose the normal resistance of the two 

junctions i?^ = 2i?„/(l — a/), = 2i?„/(l + aj) thus keeping the characteristic voltage Vc = IcRn = Iq^u — IqR^ 
and the corresponding intrawell relaxation frequency lo^) — ^Vc of each junction constant (i?„ is the SQUID normal 
resistance). This choice corresponds to an /q asymmetry due to different junction areas. We calculated the normalized 
mean voltage < v >= ^ ^ ^ V{i{ty)dt for an ac excitation i{t) = iacSa\{ijjt) with normalized amplitude iac = 
lac/Ic- The integration period T was chosen to be T > 2tt/uj (for uj < wq)- We first discuss the rocking ratchet 
case (a/ ^ 0, = 0). Here, we find < v > {iac) characteristics as shown in Fig.j3, which are strikingly similar 
to those calculated for a 1-dimensional rocked thermal ratchet For F = [Fig. |3[ (a)] we find for adiabatically 
slow excitation, uj — 0.001 (a) = lo/ujq), a lower (i;) and an upper threshold for iac- For iac < ii, no mean 
voltage appears. For i; < iac < iu the mean voltage increases monotonically with increasing iac and for iac ^ 
the voltage decreases monotonically. For oj — 0.01, the overall shape oi < v > (iac) is similar to the uj ^ case, 
however, voltage steps at < v >= ulj appear. These steps can be interpreted as Shapiro steps where the phase 
dynamics synchronizes with the external excitation. For a given iac and within one excitation period the phase " rolls" 
m2TT in one current direction and fc27r in the opposite current direction (m, k: integer). On average the phase evolves 
{m — k)2'K per excitation period and hence < v >= {m — k)u)- Within ii < iac < iu, we have k = and m increases with 
increasing iac- For iac > iu both, m and k increase with increasing iac, while (to — k) exhibits a switching behavior 
[see inset Fig. || (a)] with an overall decrease with increasing iac- As iac — > 00 we have (to — fc) ^ 0. With increasing 
frequency the value of (to — k) decreases. For cD = 0.1 a finite voltage appears in certain Zac-intervals only, where 
plateaus of equal voltage appear, i.e. (to — fc) = 0, 1. A further increase of lo confines the locking condition to obtain 
(771 — /c) > to smaller Zac-intervals. For lu =^ 0.5, without noise, no mean voltage is predicted. If one introduces small 
thermal fluctuations (F = 0.03) [see Fig. |q(b)], the condition for {m — k) > is relaxed and a stochastic resonance 
like effect sets in (for a recent review see |lq]). While the maximum voltage for each frequency is reduced the sharp 
plateaus observed for F = smear out. Furthermore, even for lu — 0.5 there exists an iac-range wherein < v >^ 0. 

We next consider the flashing ratchet case (a^ ^ 0). Here, a novel feature appears: the mean voltage undergoes 
multiple sign reversals with increasing amplitude iac as shown in Fig. IJ. The envelope of the peak voltages exhibits 
an oscillatory damped behavior, and with finite F = 0.03 the step like behavior is again smeared out. An additional 
/g-asymmetry [Fig. ^(b)] further complicates the < v > (iac)-characteristics. For F = 0, the oscillation period of the 
envelope is clearly smaller as in Fig. ^(a), however, this oscillation is less obvious if one adds small thermal fluctuations 
(F = 0.03). In any case a finite F induces further sign changes in < d > (iac)- 

To test our model we performed measurements on a YBa2Cu307 thin film dc SQUID with 24° bicrystal Joseph- 
son junctions on a SrTiOa substrate. At T =78K, we found Ic =172 ^A (F «0.04), i?„ =1.34f], Vc =232^V, 
ujo/2tt =112GHz and f3c ~ 0.07. From the estimated geometric inductance we obtain L « 30pH. However, we note 
that due to the unknown contribution from the kinetic inductance, L may be substantially larger. All measurements 
were carried out in a magnetically shielded cryostat. To apply the ac bias current we used two different techniques: 
in the adiabatic frequency regime up to 500kHz the ac current was fed via a twisted pair of wires to the SQUID. 
At higher frequency (12.31GHz) the microwave power was guided by a semi rigid coaxial cable into the cryostat. At 
the end of the coaxial cable an antenna structure was positioned about 1cm above the SQUID. In this geometry the 
microwave coupling is not well defined, which implies that both, in- and out- of phase currents are induced across the 
two junctions. The latter corresponds to a circulating current in the SQUID loop. Hence, a substantial self-field effect 
can be generated which is not present in the adiabatic experimental setup. For details of the experimental setup and 
SQUID layout see 

The measured mean voltage vs. Iac at adiabatically slow excitation (a;/27r = 50kHz) is shown in Fig. |^(a) together 
with the numerical simulation result for L =60 pH corresponding to Pl = 5.1. The measured flux shift A$ w 0.5<f>o 
(see Fig. ^ yields a junction asymmetry of aj = A/^qPl ~0.1. The good agreement is a further confirmation 
of a ratchet effect in our system. In contrast to the adiabatic case a much more complex behaviour occurs at 
lu/2tt — 12.31GHz {lu « 0.1) as shown in Fig. ||(b). The mean voltage oscillates with increasing ac amplitude and 
multiple sign reversals occur. At this frequency the ac bias is not well defined and the self-field effect can play a major 
role. In fact, the qualitative agreement with the simulation results shown in Fig. 0(b) for F = 0.03 suggests that the 
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assumption of a finite value of is reasonable. We believe that the observed behaviour can be ascribed to a flashing 
rocked ratchet system. 

In conclusion, we propose asymmetric dc SQUIDs as a model system for the study of the ratchet effect. We 
show that the voltage rectification effect in dc SQUIDs occurs due to broken reflection symmetry in the 2-dimensional 
SQUID potential and present a model which accounts for this observation by introducing a critical current asymmetry 
and a self field effect. The latter, although found to be responsible for a rich and complex dynamic behaviour, requires 
better control of high frequency coupling to the SQUID in further experiments. For the asymmetry due to different 
critical junction currents we obtain a simple rule to design SQUIDs with a large ratchet effect. It should be pointed 
out that the observed ratchet effect is no genuine property of high- Tc SQUIDs. In fact, using conventional low- Tc 
SQUIDs might have advantages due to the superior low spread junction technology. However the voltage scale on 
which the ratchet effect occurs is determined by the IcRn product which is larger for high-Tj, SQUIDs Finally 
we emphasize that the observation of a ratchet effect in asymmetric dc SQUIDs opens up the perspective for a variety 
of experimental studies in a system with straightforward detection of directed motion, simply by measuring voltage. 
The system offers experimental control over important parameters such as amplitude and frequency of the external 
force, damping coefficient and thermal noise parameter. Furthermore, the impact of different types of fluctuations 
might be interesting to investigate. 

We gratefully acknowledge R. Kleiner for providing us with his powerful simulation tool, and we thank T. Bauch 
and A. Schadschneider for fruitful discussions. 
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Fig. 1 Contour plot of the SQUID potential Ug {(3l = 1-5, fa = 0.25) with and without Iq asymmetry (upper row: 
ai = 0, lower row: aj = 0.3) and for both directions of bias current i. Arrows indicate direction of motion along 
phase trajectories (dashed hnes). 

Fig. 2 Measured dc voltage V vs. applied flux for different values of bias current / with a flux shift A$ = 
|$i — $2! ~ 0.5$0) with $1^2 defined as the flux values of corresponding extrema. 

Fig. 3 Calculated mean voltage < v > across dc SQUID with Iq asymmetry (aj = 0.3) vs. ac amplitude iac without 
thermal noise (a) and for F = 0.03 (b) for different frequencies and fa = 0.25. The inset in (a) shows an expanded 
view for a) = 0.01 at ~1- The labels for w in (b) are also valid for (a). 

Fig. 4 Calculated mean voltage < v > across dc SQUID with self-fleld effect vs. ac amplitude iac without thermal 
fluctuations (thick line) and for F = 0.03 (dots); fa = 0.25. 

Fig. 5 Mean voltage across dc SQUID vs. ac amplitude at fa = 0.25: (a) comparison of experiment and calculation 
in the adiabatic limit; (b) experiment for i2) ~ 0.1. 
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